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On Landau's prediction for large-scale fluctuation of turbulence energy dissipation 
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Kolmogorov's theory for turbulence in 1941 is based on a hypothesis that small-scale statistics are 
uniquely determined by the kinematic viscosity and the mean rate of energy dissipation. Landau 
remarked that the local rate of energy dissipation should fluctuate in space over large scales and 
hence should affect small-scale statistics. Experimentally, we confirm the significance of this large- 
scale fluctuation, which is comparable to the mean rate of energy dissipation at the typical scale 
for energy- containing eddies. The significance is independent of the Reynolds number and the 
configuration for turbulence production. With an increase of scale r above the scale of largest 
energy-containing eddies, the fluctuation becomes to have the scaling r~ 1//2 and becomes close to 
Gaussian. We also confirm that the large-scale fluctuation affects small-scale statistics. 



I. INTRODUCTION 

For locally isotropic turbulence, Kolmogoro\ii based 
his theory in 1941 on a hypothesis that small-scale statis- 
tics are uniquely determined by the kinematic viscosity v 
and the mean rate of energy dissipation (e) . Landau^ re- 
marked as follows. The local rate of energy dissipation e 
should fluctuate in space over scales of large eddies. This 
large-scale fluctuation should not be universal but should 
be different for different flows. Since the large-scale fluc- 
tuation should affect small-scale statistics, there should 
be no universal law that describes the small-scale statis- 
tics. 

Landau's remark has been wrongly regarded as a pre- 
diction for small-scale intcrmittency, i.e., enhancement 
of a small-scale fluctuation in a fraction of the volume. 
This is because, when Kolmogoro\^ revised his theory 
in 1962 to incorporate small-scale intermittency, he gave 
too much credit to LandauA^ For small scales, Landau 
remarked not about intcrmittency but about universality. 

The small-scale statistics in Landau's original remark 
are the second-order moments of velocity increments 
alone. Here we extend Landau's remark to consider the 
universality of other small-scale statistics. 

The importance of Landau's remark is not only to 
small-scale statistics. Any potentially significant phe- 
nomenon at large scales is in itself important because 
the large scales dominate the turbulence energy. 

Despite numerous studies of small-scale intermittency 
in the energy dissipationj^daiiSiL&iLiSiAi its large-scale 
fluctuation has not attracted much interest. Experimen- 
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tally, the fluctuation has been known to exist on some 
level, but it has not been discussed with any highlighting. 
The existing works are all theoretical. Oboukhov 12 con- 
sidered that the fluctuation is significant. Kraichnan 13 
considered that the fluctuation is smoothed out by spatial 
mixing, which could be caused by the action of pressure 
fluctuation^ Frisch^ considered that the fluctuation is 
significant in some specific configurations for turbulence 
production. Thus, even for the significance of the fluctu- 
ation, no consensus has been reached on. 

We experimentally assess Landau's remark, i.e., 
whether the large-scale fluctuation of energy dissipation 
is significant, whether the large-scale fluctuation is uni- 
versal, and whether the large-scale fluctuation affects 
small-scale statistics. The assessment requires more than 
one type of turbulence. We study grid and boundary- 
layer turbulence at several Reynolds numbers. With a 
hot-wire anemometer and Taylor's frozen-eddy hypothe- 
sis, we obtain a one-dimensional cut of the velocity field. 
The local rate of energy dissipation per unit mass is ob- 
tained along the streamwise position x as 



du(x) 
dx 



I5v 



dv(x) 



dx 



(1) 



Here u and v are the streamwise and transverse ve- 
locities. These rates are surrogates of the true rate 
v ^2i j(di u j + 9jUi) 2 /2, where i and j denote coordinate 
axes. Their averages are nevertheless exact if turbulence 
is locally isotropic. The energy dissipation at the scale r 
is obtained by coarse-graining the local rate^ 



s(r,x) 



1 



-r/2 



e(x + x')dx' . 



(2) 



-r/2 



For this and other equations, the subscript u or v is omit- 
ted if it is not necessary. The fluctuation of the energy 
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dissipation e(r,x) around its average (e(r, x)) — (e(x)) is 
studied as a function of scale r. 

We study the features expected for largest scales in Sec. 
ITT1 Grid turbulence is studied in Sec. IIIII Boundary-layer 
turbulence is studied in Sec. IIVI The dependence on 
the Reynolds number is studied in Sec. [Vj We summa- 
rize our conclusions and discuss the effect on small-scale 
statistics in Sec. I VII Taylor's frozen-eddy hypothesis for 
large scales is studied in Appendix. 



II. LARGEST-SCALE FEATURES 

Suppose that the local rate of energy dissipation e(x) is 
obtained on a one-dimensional cut of a turbulence flow. If 
the turbulence is homogeneous along the one-dimensional 
cut, the correlation function e (r) at the scale r is 

Mr) = Mx + r) - (e(x))}[e(x) - (e(x))}). (3) 

Here (■} denotes a spatial average over the position x. 
The correlation length L e is 



L F = 



f °° Mr)dr 
MV ' 



(4) 



We assume that the correlation is insignificant or absent 
above a certain scale r*, i.e., <j> e (r)/(j> E (0) ~ at r > r* > 
L e . This assumption is natural for the correlation of any 
quantity in real turbulence that only has a finite extent. 

The correlation function <f> £ (r) allows us to obtain the 
mean-square fluctuation of the energy dissipation e(r, x) 
around its average (e(r,x)): 

a £ {rf = ([e(r,x) - (e(r,x))} 2 ) = 1 f \r -r')Mr')dr' '. 



If the scale r is much greater than the scale r* , we have 

2L £ M0) 



cr e (r 



(6) 



Thus the root-mean-square fluctuation scales as o e (r) oc 

r -l/2_ 

We explain the scaling u e {r) oc r -1 / 2 . If the scale r 
is much greater than the scale r*, the energy dissipation 
e(r, x) is written with the sum over subregions that have 
a width r*: 



r/r* 



e(r,x) 



(7) 



where x n = x — (r + r*)/2 + nr t . Since the energy dis- 
sipations e(r*,x n ) are mutually independent, the mean- 
square fluctuation for J2 n £ { r *i x n) is J2n a s( r *) 2 - Then 
we have 



Mr) 



2 T# ®e (y* ) 



(8) 



which again leads to the scaling a E (r) oc r^ 1 / 2 . Thus this 
scaling is not dynamical but statistical. 

There is another largest-scale statistical feature. The 
total number of subregions r/r* in Eq. Q is large. Then 
the central limit theorem ensures that the energy dissipa- 
tion s(r,x) obeys a Gaussian distribution at least in the 
vicinity of the average (e(r, x))r& provided that the fluc- 
tuation of energy dissipation a e (r) is sufficiently smaller 
than the average (e(r,x)}. 

The scaling u e {r) oc r -1 / 2 and the tendency for Gaus- 
sianity imply that rs(r, x) is additive at largest scales. 
It is not additive at the smaller scales. The value of an 
additive quantity for a region is the sum of its values 
for the subregions that are mutually independent [Eq. 
1(7)1]. Many examples exist in the statistical mechanics 
and thermodynamics^ There the scales of interest are 
much greater than the scales where the correlation is sig- 
nificant. 



III. GRID TURBULENCE 
A. Experiment 

The experiment was done in a wind tunnel of the Me- 
teorological Research Institute. We use the coordinates 
x, y, z in the streamwise, spanwise, and floor-normal di- 
rections. The corresponding turbulence velocities are u, 
v, and w. The origin i = y = z = 0m is taken on the 
tunnel floor at the entrance to the test section. Its size 
was Sx = 18 m, Sy = 3 m, and Sz — 2 m. 

We placed a grid across the entrance to the test section 
at x = m. The grid consisted of two layers of uniformly 
spaced rods, the axes of which were perpendicular to each 
other. The separation of the axes of adjacent rods was 
0.20 m. The cross section of the rods was 0.04 x 0.04 m. 
We set the mean wind to be U ~ 10ms -1 . 

The streamwise (u) and spanwise (v) velocities were 
simultaneously measured using a hot-wire anemometer. 
The anemometer was composed of a constant tempera- 
ture system and a crossed-wire probe. The wires were 
made of platinum-plated tungsten, 5 /xm in diameter, 
1.25 mm in sensitive length, 1mm in separation, and 
280 °C in temperature. 

The measurement was done on the tunnel centerline 
at x = 4m, where the flatness factors (u(x) 4 ) / (u(x) 2 ) 2 
and (v(x) 4 )/ (v(x) 2 ) 2 were close to the Gaussian value of 
3. Thus turbulence was well developed, so eddies with 
various sizes and strengths filled the space randomly and 
independently^*^ If the measurement position had been 
too close to or far from the grid, turbulence should have 
been still developing or already decaying, and its flat- 
ness factors should have been different from the Gaussian 
value. The ratio (u 2 ) 1 ^ 2 / (v 2 ) 1 ^ 2 was close to unity. 

The signal was linearized, low-pass filtered at 15 kHz 
with 24 dB per octave, and then sampled digitally at 30 
kHz with 16-bit resolution. The total length of the signal 
was 3 x 10 s points. 
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TABLE I: Summary of flow parameters. 



Quantity Units Grid Boundary layer 











z — 0.25 m 


z — 0.70 m 


JVlcan strcamwisc velocity 


U 


ms -i 


10 57 


7 05 


10 22 


Strc 3,m wise velocity fi.uctu.tit ion 


(u 2 ) 1/2 


m8 _! 


0.525 


1.48 


0.642 


Spairwisc velocity fluctuation 


(v 2 ) 1 ' 2 
\ u 1 




O.515 


1.22 


0.544 


Strcamwisc 11 fit noss factor 






3.02 


2.71 


7.83 


Spanwise flatness factor 


\ u II \" 1 




2.99 


3.02 


7.36 


Air temperature 




°C 


11.6-12.0 


14.9-15.9 


14.4-14.9 


Kinematic viscosity 


// 


cm 2 s — 1 


0.142 


0.145 


0.145 


Mean energy dissipation rate {c u ) 


(e„) = ISvMd^u) 2 ) 


m 2 s- 3 


1.22 


5.45 


0.575 


Mean energy dissipation rate (e v ) 


<£„> = i5^<(a X D) 2 >/2 


2 —3 

m s 


1.20 


4.05 


0.493 


Correlation length (u) 


L„ = / o -0„(r)dr/0 u (O) 


cm 


17.4 


42.3 


35.9 


Correlation length (v) 


L v =/~0„(r)dr/^(O) 


cm 


4.43 


5.78 


7.60 


Correlation length (e u ) 


Leu = Jo°° ' <t>e*(r)dr/fa u (0) 


cm 


0.545 


0.970 


2.69 


Correlation length {e v ) 


L ev = /« <fe„(r)<ir/0 e „(O) 


cm 


0.414 


0.753 


2.03 


Taylor microscalc 


\=[2{v 2 )/({d !c v) 2 )} 1 ' 2 


cm 


0.688 


0.896 


1.14 


Kolmogorov length 


I) = (, 3 /<^» 1/4 


cm 


0.0221 


0.0166 


0.0280 


Microscalc Reynolds number 






249 


756 


428 



In addition, to study the average and fluctuation of 
energy injection at x — 4 m, the streamwise velocity was 
simultaneously measured at x = 3.75 and 4.25 m. We 
used two single-wire probes. The total length of the sig- 
nal was 10 8 points. The other conditions were the same 
as in the above measurement. 

The flow parameters at x = 4 m are summarized in Ta- 
bleHJ We used Taylor's frozen-eddy hypothesis to convert 
temporal variations into spatial variations in the stream- 
wise direction (see Appendix). The velocity derivatives 
were obtained from finite differences, e.g., 



du[x) u(x + Sx) — u[x — 5x) 



dx 



28x 



(9) 



Here 5x is the sampling interval. Since the sampling 
interval was small enough, the higher-order accuracy for 
the velocity derivatives is not necessary.^ 



B. Results and discussion 

The statistics of energy dissipation and relevant quan- 
tities are studied over a wide range of scale r. Through- 
out the r range, statistical significance is satisfactory be- 
cause our data are long. 

Figure^a) shows the correlation functions for the local 
rates of energy dissipation e u (x) and e v (x) and for the 
velocities u(x) and v(x). The correlation lengths are also 
shown. For the velocities, we defined as 



4> u (r) = (u(x + r)u(x)) and L u = - " , 

Mr) = (v(x + r)v(x)) and L v =- ^° ^ v ^ dr 



(10a) 



0,(0) 



. (10b) 



These correlation functions and correlation lengths of- 
fer information about the scale structure of turbulence. 



The velocity correlations 4> u {f) and 4> v {r) exist up to the 
scale of largest energy-containing eddies. The correla- 
tion length L u corresponds to the mean scale of energy- 
containing eddies. Since the local dissipation rates be- 
long to small scales, their correlation functions <p eu (r) 
and (j) Ev (r) decay fast (see Ref. for another interpreta- 
tion) . They nevertheless exist up to about the correlation 
length L U &H 

Figure ^b) shows the root-mean-square fluctuations 
of energy dissipation cr Su (r) and a Sv (r). They are en- 
hanced at small scales owing to intermittency. They are 
still comparable to the mean energy dissipation (e(x)} 
at about the correlation length L u . Thus, we confirm 
Landau's remark^ that the energy dissipation should sig- 
nificantly fluctuate over scales of energy-containing ed- 
dies. Consistent results are seen, albeit not explicitly 
stated, in past experimental worksiS-ii The statistical 
scaling a e (r) oc r -1 / 2 starts at about the scale where 
the velocity correlations <p u {f) and 4>v{r) vanish, i.e., the 
scale of largest energy-containing eddies. 

Figure^b) also shows the flatness factor for the fluctu- 
ation of the energy dissipation s(r,x) around its average 
(e{r,x)): 



flatness factor 



([e(r,x)-(e(r,x))]*) 
{[e(r,x) - (e(r,x)}} 2 ) 2 



(11) 



The flatness factor is enhanced at small scales owing to 
intermittency. With an increase of scale above the cor- 
relation length L Ul the flatness factor approaches the 
Gaussian value of 3. The exact value is not achieved 
even at largest scales^ although the discrepancy is too 
small to be discernible in our diagram. This is because 
the central limit theorem does not necessarily lead to an 
exactly Gaussian distribution. The discrepancy from the 
Gaussian value would be larger in the higher-order statis- 
tics, the reliable computation of which requires the longer 
data. 
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10° 10 1 10 2 10 3 10 4 10 5 10 6 
scale (r) 
Kolmogorov length 



FIG. 1: Statistics for grid turbulence as a function of 
scale r normalized by the Kolmogorov length, (a) Corre- 
lation functions 0e„(r)/^ Su (O), <t>s v (r)/(j)e v (0), <j>u(r) / t/> u (0) , 
and 4>v(r) / <j> v (0). The arrows indicate the correlation lengths 
L Su , L Bv , L u , and L v . (b) Root-mean-square fluctuations 
a> u (r) and o Ev (r). They are respectively normalized by 
(e u {x)) and (e v (x)). The dotted line indicates the r _1//2 scal- 
ing. We also show flatness factors for e u (r,x) and e v (r,x). 
The dotted line indicates the Gaussian value of 3. (c) Av- 
erages and root-mean-square fluctuations of energy injection, 
transfer, and dissipation. They are respectively normalized 
by (e„(a;)), {e v {x)), and (e u (x)). The dotted lines indicate 
the r -1 / 2 and r~ 3 ^ 2 scalings. 

Except at smallest scales, the statistics for the sur- 
rogate dissipation rates e u (x) and e v {x) are almost the 
sameiii It is thereby expected that the true dissipation 
rate would yield almost the same statistics. 

Figure Q^c) shows quantities in the scale-by-scale en- 
ergy budget equation for decaying homogeneous isotropic 
turbulence >iLi& 

5{Su(r,x) 3 ) Vbv d{Su(r : x) 2 ) 



Here 8u(r, x) is the velocity increment u{x + r) — 
u{x). We coarse-grained the local energy injection 
— 15 J Q dt(Su 2 )r' dr' /4r 5 , the transfer — 5<5u 3 /4r, and the 
dissipation \hvd r {5u 2 ) /2r over the scale r and computed 
their averages and root-mean-square fluctuations around 
the averages (for the energy injection, see also Appendix). 



The mean energy transfer is close to the mean energy dis- 
sipation (s(x)} in the inertial range, which extends up to 
the mean scale of energy-containing eddies. There the 
fluctuations of energy injection and transfer are compa- 
rable to each other and are several times greater than 
the mean energy dissipation. In particular, the fluctua- 
tion of energy injection is maximal. These fluctuations 
become to have the statistical scalings r -1 / 2 and r -3 / 2 , 
respectively, 19 at the scale where the fluctuation of en- 
ergy dissipation cr e (r) becomes to have the statistical 
scaling r -1 / 2 . 

The energy dissipation occurs at the end of the en- 
ergy cascade and hence depends on the scale-by-scale 
energy transfer. 13 Although the mean energy transfer is 
downward, the local energy transfer is as often upward as 
downward at each scale because its fluctuation is strong 
[Fig. (He)]. This fluctuation of energy transfer causes the 
large-scale fluctuation of energy dissipation. Since the 
former is stronger than the latter [Figs, ^c) and^b)], 
spatial mixing is at workji 3 - albeit not completed The 
fluctuation of energy transfer is associated with the fluc- 
tuation of energy injection. Since the local energy injec- 
tion is as often negative as positive [Fig. ^c)], the indi- 
vidual energy-containing eddies as often gain as lose their 
energies. There is accordingly the significant fluctuation 
of energy dissipation over scales of energy-containing ed- 
dies [Fig. IHb)]. If the scale r exceeds the scale of largest 
energy-containing eddies, the fluctuations of energy in- 
jection, transfer, and hence dissipation reflect the ran- 
dom and independent distribution of energy-containing 
eddies. They have the statistical scalings r -1 / 2 , r~ 3 / 2 , 
and r~ x l 2 [Figs, ^c) and^b)]. The fluctuation of en- 
ergy dissipation is close to Gaussian [Fig. IHb)]. 



IV. BOUNDARY-LAYER TURBULENCE 
A. Experiment 

The experiment was done in the same wind tunnel with 
the same instruments as for the grid turbulence. We de- 
scribe specific features. The flow parameters are summa- 
rized in Tabled 

Over the entire floor of the test section of the wind 
tunnel, we placed blocks as roughness for the boundary 
layer. The block size was 8x = 0.06 m, Sy = 0.21m, and 
Sz = 0.11m. The spacing of adjacent blocks was 5x = 
Sy = 0.5 m. We set the incoming wind to be 10ms -1 . 

The measurement positions were at x = 12.5 m, where 
the boundary layer was well developed. The 99% thick- 
ness, i.e., the height at which the mean streamwise veloc- 
ity U is 99% of its maximum value U, was 0.77 m. The 
displacement thickness f Q (1 — U/U)dz was 0.20 m. Here 
z is the height for the velocity U ^ The log-law sublayer 
was at z ~ 0.15-0. 35m. 

The measurement was done at z = 0.25 and 0.70 m. 
The height z = 0.25 m was in the log-law sublayer, 
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FIG. 2: Same as in Fig. 0but for boundary-layer turbulence 
at z = 0.25 m. 
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FIG. 3: Same as in Fig. 0but for boundary-layer turbulence 
at z = 0.70 m. 



where various eddies filled the space randomly and 
independently^ The flatness factor (v(x) 4 } / (v(x) 2 ) 2 was 
close to the Gaussian value of 3. The height z = 0.70 m 
was in the wake sublayer, where eddies did not fill the 
space. They intermittently passed the probe. 

Since turbulence was not weak compared with the 
mean streamwise velocity, the u component at smallest 
scales measured by the crossed-wire probe was contami- 
nated with the w component that is perpendicular to the 
two wires of the probe. The v component was free from 
such contamination. 21 Hence we prefer the rate e v (x) to 
the rate e u (x). 

Turbulence was not isotropic at large scales, but it 
should have been isotropic at small scales. This is be- 
cause, with a decrease of scale r, the observed ratio be- 
tween 2r{Sv(r,x) 2 ) and d r (r 2 {5u(r,x) 2 )) once becomes 
the isotropic value of unity at about the correlation 
length L u . 



B. Results and discussion 

Figures [21 and [21 show the same statistics as Fig. ^ 
except for the energy injection that is due to large-scale 
shear of the boundary layer. There is no known formula 
to compute the local injection of energy budget in such a 
shear flowiSS The energy transfer and dissipation are still 
of physical significance and their averages satisfy Eq. 1|12|) 
at small scales where the energy injection is negligible and 



the turbulence is isotropic (see Sec. IIV All 

The mean energy transfer is different. At z = 0.25 m 
[Fig. EJc)], the inertial range is wide. At z = 0.70m 
[Fig. Etc)], the inertial range is narrow. There were 
only few eddies. Their main bodies were at z <^ 0.70 m. 
Their small streamwise sizes at z — 0.70 m correspond 
to the upper edge of the inertial range. The correlation 
length L u is relatively large because it reflects the spatial 
distribution of those eddies. 

Despite the difference of the mean energy transfer in 
Figs. [21 and [21 we find common features, albeit not ex- 
actly universal. They were also found in Fig. ^for grid 
turbulence. The local rate of energy dissipation has a 
correlation up to about the correlation length L u . There 
the fluctuation of energy dissipation a e (r) is comparable 
to the mean energy dissipation (e(x)). The fluctuation of 
energy transfer is stronger. These fluctuations become to 
have the statistical scalings r -1 / 2 and r~ 3 / 2 at about the 
scale where the velocity correlations vanish. The fluctu- 
ation of energy dissipation becomes close to Gaussian. 

Frischi*^ considered that the large-scale fluctuation of 
energy dissipation is significant if flow parameters change 
in space over a scale larger than the mean scale of energy- 
containing eddies. An example is that turbulence does 
not fill the spaced This is the case at z = 0.70 m. How- 
ever, not only at Z = 0.70 m [Fig. |2tb)] but also at 
z = 0.25 m [Fig. Hfb)], the large-scale fluctuation of en- 
ergy dissipation is significant. 
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FIG. 4: Root-mean-square fluctuation a Cv (r) as a function of 
scale r. The ordinate is normalized by (e v (x)). The abscissa 
is normalized by the correlation length L u . The dotted line 
indicates the r -1 / 2 scaling, (a) Grid turbulence. The data 
for ReA = 249 are from our present experiment. The data 
for Re A = 105, 165, 225, 292, and 329 are from our past 
experiment!— (b) Boundary layer. The data for Re A — 756 
are from our present experiment. The data for ReA = 295, 
430, 655, 861, 1054, and 1258 are from our past experiment. 24 
The inset shows the dependence on the microscale Reynolds 
number Re a at r = L u . 



V. DEPENDENCE ON THE REYNOLDS 
NUMBER 

The dependence of the fluctuation of energy dissipation 
<j Sv (r) on the microscale Reynolds number ReA is studied 
in Fig. 21 We use the data from our present experiments 
and also those from our past experiments^^ 

Figure 0Ja) shows the fluctuation <r Sv {r) for grid tur- 
bulence at Re A = 105-329. The grid or the mean 
wind velocity was not the same, but the flatness factors 
{u{x) A ) I (u(x) 2 ) 2 and (v(x) A ) / (v(x) 2 ) 2 were close to the 
Gaussian value of 3. The ratio (u 2 ) 1 / 2 / (v 2 ) 1 / 2 was close 
to unityi^ 

Figure^Jb) shows the fluctuation a £v (r) for boundary- 
layer turbulence at ReA = 295-1258. The data were ob- 
tained in log-law sublayers at the same streamwise po- 
sition x over the same roughness for different incoming- 
wind velocities. The flatness factor (v(x) A ) / (v(x) 2 ) 2 was 
close to the Gaussian value of 3iSi 



Regardless of the Reynolds number, the fluctuation 
a Bv (r) at the correlation length L u is comparable to the 
mean dissipation rate (e v (x)). The boundary-layer tur- 
bulence has a trend that the fluctuation a £v (r) is smaller 
at a higher Reynolds number. However, at ReA = 9000 in 
an atmospheric boundary layer the fluctuation a Sn (r) 
at the correlation length L u is still comparable to the 
mean dissipation rate (e u (x)). The trend at ReA 3> 1000 
is insignificant or absent. We could attribute the trend to 
spatial mixing during the energy cascade>ii2^ The cas- 
cade depth and duration are significantly greater at a 
higher Reynolds number so far as it is not very high. For 
the grid turbulence, the trend is not discernible because 
the Reynolds number does not span a wide range. 



VI. CONCLUDING DISCUSSION 

Using grid and boundary-layer turbulence, for the first 
time, we experimentally confirm Landau's remark 2 that 
the energy dissipation should significantly fluctuate in 
space over large scales. This large-scale fluctuation is 
caused by the large-scale fluctuation of scale-by-scale en- 
ergy transfer. Although the large-scale fluctuation of 
energy dissipation is smoothed by spatial mixing^ the 
smoothing is not completed 

The large-scale fluctuation of energy dissipation is not 
exactly universal as remarked by Landau^ Nevertheless, 
there are features that are independent of the Reynolds 
number and the configuration for turbulence production. 
The fluctuation a £ (r) is comparable to the mean energy 
dissipation (e(x)} at the correlation length of streamwise 
velocity L u . If the scale r exceeds the scale where the 
velocity correlations vanish, the fluctuation has the sta- 
tistical scaling a E (r) cx r -1 / 2 and is close to Gaussian. 

The correlation length L u corresponds to the mean 
scale of energy-containing eddies. The velocity correla- 
tions exist up to the scale of largest energy-containing 
eddies. Hence, although the large-scale fluctuation of en- 
ergy dissipation is caused by that of energy transfer, they 
are ultimately caused by individual energy-containing ed- 
dies. Exceptionally, when turbulence does not fill the 
space, the large-scale fluctuations are determined by the 
distribution of energy-containing eddies. 

There are global quantities that exhibit significant 
temporal fluctuations. An example is the total energy 
injection required to sustain the Karman flow, i.e., tur- 
bulence between two counter-rotating disks— i*2L2& The 
fluctuation is considered to be associated with the fluctu- 
ation of energy transfer. When the system size, e.g., the 
diameter of the two disks relative to the separation be- 
tween them, is much greater than the correlation length, 
the fluctuation is Gaussian. The existence of such fluctu- 
ations is consistent with our results for large-scale spatial 
fluctuations. 

Landau^, also remarked that the large-scale fluctuation 
of energy dissipation should affect small-scale statistics 
in the inertial and dissipation ranges. We finally discuss 
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103 
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FIG. 5: Exponent 7 for the inertial-range scaling a e „ (r) oc 
r ' as a function of the microscale Reynolds number ReA- 
The filled circles denote grid turbulence. The open circles 
denote boundary-layer turbulence. The data are the same as 
in Fig. [1 



and hence the large-scale fluctuations are different. The 
exponent sequence is accordingly different. Consistent 
results exist for the scaling tf> eu (r) oc and the depen- 
dence of 5u(r,x) and e(r, x) on a large-scale quantity, 
M ^ | 30 i 33 4 34 J 35 Thus, since the large-scale fluctuations 
are not exactly universal, the high-order small-scale 
statistics are not universal, although some conditional 
statistics might be universalii2i2£i2£ 

Having confirmed Landau's remark, we underline that 
turbulence is much more fluctuating over large scales 
than it is usually assumed to be. The individual energy- 
containing eddies as often gain as lose their energies. The 
local energy transfer is as often upward as downward. 
These significant fluctuations lead to the large-scale fluc- 
tuation of energy dissipation and also affect small-scale 
statistics. 



this remark. 

The large-scale fluctuation of energy dissipation and 
that of energy transfer do not affect the small-scale statis- 
tics considered by Kolmogorov^ i.e., the mean energies 
(6u(r,x) 2 ) and (Sv(r, x) 2 ). They are determined by the 
averages of energy transfer and dissipation, which satisfy 
the energy budget equation (|12|) in spite of their strong 
fluctuations [Fig. ^c)-|3[c)]. The Kolmogorov constant 
(5u(r, x) 2 )/(r{e(x))) 2 / 3 for the inertial range is actually 
universal among various flows' When turbulence does 
not fill the space, the Kolmogorov constant is different, 
but the universal value is obtained if we consider the tur- 
bulence regions alone^ In addition, the energy spectrum 
for small scales is universal if it is normalized with the 
kinematic viscosity v and the mean energy dissipation 

The large-scale fluctuation of energy dissipation and 
that of energy transfer affect high-order small-scale 
statistics, e.g., (5u(r,x) n ) for n > 4 and (s(r, x) n ) for 
n > 2. They represent the fluctuations of energy, en- 
ergy transfer, and energy dissipation, which suffer from 
the large-scale fluctuations. For example, (Su(r, x) 6 ) cor- 
responds to the mean-square fluctuation of the energy 
transfer — 5<5w 3 /4r. The effect of the large-scale fluctua- 
tions is not always significant. We are still able to demon- 
strate its existence using an inertial-range scaling in Figs. 

una 

cr ev (r) oc r~ 7 , (13) 

which is more significant than the famil- 
iar scalinga 3 ! 4 . 5 . 6 . 7 . 8 . 9 . 10 ! 11 ! 18 (Su(r,x) n ) oc r^, 
(e(r,x) n ) oc r T " , and 4> e (r) oc r~ M at least in our 
data (see also Ref. 0). The data for Fig. 0] are used 
to compute the exponent 7 over the r range from 30?7 
to 0.3L„. Figure [3] shows the results as a function 
of the microscale Reynolds number ReA- The grid 
turbulence and the boundary-layer turbulence make 
up different sequences^ Although the flatness factor 
(v(x) 4 ) I (v(x) 2 ) 2 was commonly close to the Gaussian 
value of 3, the configuration for turbulence production 
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APPENDIX: TAYLOR'S HYPOTHESIS 

Taylor's frozen-eddy hypothesis leads to the conversion 
of the streamwise position x L and time t L in the labora- 
tory frame into the streamwise position x T and time t T 
in a virtual reference frame (hereafter, the Taylor frame): 

x T = -Ut L and t T = (A.l) 

Here U is the mean streamwise velocity. The second con- 
version is not familiar but is necessary because eddies are 
not exactly frozen (see below). Then a stationary signal 
is converted into a streamwisely homogeneous signal. A 
streamwisely homogeneous signal is converted into a sta- 
tionary signal. 

Taylor's hypothesis is valid provided that the turbu- 
lence strength (u 2 ) 1 / 2 /U is small enough. This is the case 
in our experiment (Table fU- For turbulence strengths 
close to ours, the validity of Taylor's hypothesis was 
confirmed using data obtained simultaneously with two 
probes separated by streamwise distances^ 

The Taylor frame is locally identical to a reference 
frame moving with the mean stream. With an increase of 
scale, Taylor's hypothesis becomes a mere convention to 
describe measured temporal variations in terms of spa- 
tial variations. The signal in the Taylor frame is still 
expected to be consistent with some realistic turbulence. 
We thereby applied Taylor's hypothesis to all the scales. 
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Since our laboratory flows were stationary, the signals 
in the Taylor frame represent flows that are homogeneous 
in the streamwise direction. We actually obtained the 

— 1/2 

scaling a E (r T ) oc r T and the tendency for Gaussianity 
expected for largest scales in homogeneous turbulence. 

Since our laboratory flows were inhomogeneous in the 
streamwise direction, the signals in the Taylor frame rep- 
resent nonstationary flows. We demonstrate this fact for 
the grid turbulence. If turbulence is isotropic, homoge- 
neous, and decaying, the energy budget equation isiLiSi 

15 f r T d(6u 2 ) , 4 , 5(5u 3 ) lbvd{5u 2 ) _ / x 

"44 J ~dt^ rr T ~^r + 2^^~~ {£} - 

(A.2) 

This is identical to Eq. (|12[l . The first term in the left- 



hand side is from the nonstationarity and represents the 
injection of energy budget. We estimated it from the 
data at x L — 3.75 and 4.25 m. The other terms were 
estimated from the data at x L = 4 m. They are com- 
pared in Fig. dc). The energy-budget equation l|A.2|) is 
satisfied throughout the scales. 

The above estimation appears to imply that the en- 
ergy budget is supplied from the inhomogeneity of tur- 
bulence along the streamwise direction in the laboratory 
frame. However, once we have used Taylor's hypothe- 
sis, it is important to distinguish the laboratory frame 
from the Taylor frame. Grid turbulence in the Taylor 
frame is homogeneous in the streamwise direction and 
nonstationary. This nonstationarity supplies the energy 
budget. 
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